Quantum dots in magnetic fields: thermal response of broken symmetry phases 
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We investigate the thermal properties of circular semiconductor quantum dots in high magnetic 
fields using finite temperature Hartree-Fock techniques. We demonstrate that for a given magnetic 
field strength quantum dots undergo various shape phase transitions as a function of temperature, 
and we outline possible observable consequences. 



t-H . The quantum mechanical characteristics (i.e. ground- and excited-state energy levels) of semiconductor quantum 
dots were recently investigated as a function of increasing magnetic field and electron number 0J^] . Unique structures 
that break rotational symmetry within a quantum dot and develop as a function of increasing magnetic field were also 
reported in theoretical investigations (^] and associated with experimental observations of charge-density redistribution 
In this paper, we explore the thermal characteristics of the broken-symmetry phase of the charge-density 
+3 \ distribution, showing that temperature can be used to induce transitions between phases of differing spatial symmetry. 

Thermal properties of, and transitions in, strongly correlated quantum mechanical systems have been studied in 
various contexts. Perhaps the most obvious are the normal-conducting to super-conducting transitions in condensed 
matter systems. Many nuclear systems exhibit an intrinsic deformation in their ground states, and deformed-to- 
i 1 spherical transitions have been studied in many models (see for example Refs. As a result, observable conse- 

*^ . quences, such as enhanced moments of inertia, decreasing quadrupole transition moments, and decreasing correlated 
pair-transfer amplitudes may be investigated. 

In studies of many-body phenomena in quantum dots, experimental efforts have focused on mapping the magnetic 
field dependence of their ground-state structure by measuring the chemical potential via capacitance spectroscopy 
H . Cusps and steps in the chemical potential were found to clearly separate different ranges of magnetic fields HQ . 
These features were identified with phase transitions in the charge density of the quantum dot. At magnetic field 
strengths on the order of a few tesla, all electrons become spin-polarized initiating the maximum density droplet 
(MDD) phase ||, in which the density is constant and homogeneous at the maximum value that can be reached in 
s ! ■ the lowest Landau level. The stability of the MDD is determined by a competition among the kinetic energy, external 
confinement, the Coulomb repulsion between electrons, and the attraction created by the Coulomb exchange term. 
For increasing magnetic field, the charge-density distribution of the droplet reconstructs B w ith a ring of electrons 
breaking off from the MDD phase. This edge reconstruction has been shown via mean-field [|10| and density functional 
theory j| calculations to result from a rotational symmetry-breaking phase transition from the MDD to a Wigner 
molecule or Wigner crystal phase. These calculations are in good qualitative agreement with recent experimental 
results where instabilities of the MDD state and other transitions in the high magnetic field region were accompanied 
by a redistribution of the charge density [Q . 

Most of the experimental investigation of many-body phenomena in quantum dots was performed at very low 
temperatures on the order of 100 mK and focused on the evolution of the ground state of the quantum dot as a 
function of increasing magnetic field and electron number JH],[^] . The thermal characteristics of these systems are 
largely unexplored, especially for large values of the magnetic field. We explore the thermal response of the broken- 
• i-h , symmetry phase at constant magnetic field and electron number and demonstrate that transitions among a number 
' of distinct, broken-symmetry phases occur as the temperature is increased. 

In describing the ground-state and low-lying (intra-band) excitations of the iV-electron semiconductor nanostruc- 
tures, it is often sufficient to restrict consideration to the conduction band using the effective-mass approximation 
p3[ . We consider the problem of N electrons of effective mass m* in a plane, (x, y), confined by an external parabolic 
potential, V(r) = ^m*ujf ) r 2 , and subject to a strong magnetic field B = B$e z . We consider the Zeeman splitting but 
neglect the spin-orbit interaction. The Hamiltonian for such a system is 
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where the vector potential is A(fi) — {Bq/2){— y^ Xi, 0), g* — 0.54, e = 12.9, m* = 0.067m e , and huj — 3 meV. 

We solve this equation at the finite-temperature Hartree-Fock level. This is the initial step to account for electron- 
electron correlations within the quantum dot and suffices for a qualitative discussion of phase transitions and thermal 
behavior. Hartree-Fock is implemented across nuclear and atomic physics as a first step towards the solution of 
the quantum many-body problem flifl . Since the Hartree-Fock equations are nonlinear, the self-consistent potential 
obtained from their solution will not necessarily show the same symmetries as the original Hamiltonian. Of course, 
the exact wave will maintain the symmetries of the original Hamiltonian. As a point of reference, an estimate of the 
relative uncertainty in total angular momentum at high magnetic field strengths was shown to be roughly 10% by the 
authors of Ref . |L(J • 

The equations describing the static mean field (i.e., Hartree-Fock) at finite temperature for an iV-electron system 
are 

[K + W([3)-E a ]\^ a (l3)=0 (2) 

where E a are the single particle energies associated with the single-particle wave functions 4> a , and K is the one-body 
kinetic-energy operator. The one-body Hartree-Fock field, W(f3), is obtained in terms of the Hamiltonian (]l|) and 
the one-body density matrix, p(fi) as, 

W{(3) = Tr[p(f3)H] . (3) 

The single-particle wave functions are denoted as \<f> a (/?)), and (3 — 1/kT is the inverse temperature. The representa- 
tion of the one-body density is given in terms of the occupation numbers n a (f3), p(/3) = J2 a n a {[3) \ 4> a ((3)}((j) a ([3) \, 
and the Fermi-Dirac occupation numbers at chemical potential p are n a ((3) = 1/ \e^^ Ea ^^ + l] . Note that Tamura 
and Ueda jl5| studied the number fluctuations ((SN) 2 ) of a quantum dot as a function of N for different field strengths 
and found that the fluctuations were rather small (on the order of 0.1) in finite-temperature Hartree-Fock applica- 
tions. We expect qualitative agreement with experiment as Hartree-Fock is known to overestimate the magnetic field 



strength at which transitions occur |16 1JJ3J. However, for a fixed field strength, experimental evidence of the thermal 
transitions should still be qualitatively visible. 

We use a Fock-Darwin basis expansion to solve the finite-temperature Hartree-Fock equations. Since we use a high 
(« 12 T) magnetic field, we consider only angular momentum states with the n = principal quantum number. The 
electrons carry spin, and so our states are labeled by k = {Ik, Sk}, where Ik is the angular momentum projection of 
the k-th state and Sk is the spin of that state. We found convergence using fifty states for the N < 8 systems. We 
also checked our zero-temperature results with other publications fiofl for various numbers of electrons in the dot and 
found satisfactory agreement. 

We begin the discussion of our results by investigating the electron charge, angular momentum, and spin densities 
as a function of increasing temperature for the N = 6 system at Bq = 12.15 T. We show densities at representative 
temperatures of 3.87 K (the low temperature limit, v = 6), 11.97 K (before the first phase transition, v = 6), 13.65 K 
(in the second phase, v = 5), and 14.32 K (in the third phase, v = 4), where v is the number of definable high-density 
regions (or vortices) in the charge density plots. The density, shown in Fig. |l|a-d, begins as a fairly well-defined 
Wigner crystal at 3.87 K, which exhibits some degree of thermal broadening at 11.97 K. The v — 5 and v = 4 phases 
continue to show a similar amount of density in the remaining vortices, while the density of the thermally dissipated 
vortices have effectively spread through the entire dot. The angular momentum along the i?-field direction, shown in 
Fig. ^3-h, exhibits well-defined structures at low temperatures which tend to decrease rapidly as one moves through 
the various phases. Although the high charge density regions in the v = 4 phase are still well defined, the angular 
momentum in this high-temperature phase has nearly washed out. Finally, we show in Fig. [l]i-l the spin density 
defined as p s (x,y) = [p]{x, y) — pi(x,y)] / [p^{x,y) + Pi{x,y)], where p| {p\) refer to the spin up (down) density. At 
these temperatures, little appreciable spin depolarization occurs and the spin density remains above 0.8 for the entire 
region where there is appreciable charge density. 

The suddenness of the phase transitions seen in Fig. [l] become quite evident when the internal energy of the quantum 
dot is plotted as a function of the temperature. We show the three phases of the dot in Fig. ||a. Note that the v = 6 
phase exists as an excited configuration when the most probable Hartree-Fock solution is the v = 5 phase. Similarly, 
the v = 4 configuration exists as a possible excited configuration of the system even at fairly low temperatures. The 
specific heat, C v = d(H)/dT (with T in units of eV), is shown in Fig. ||b. Clearly, when the energy undergoes a phase 



2 



transition, the specific heat shows a sharp structure. This occurs since the energy is piece-wise continuous along the 
three phases. 

These calculations suggest that the quantum dot exhibits a band structure of many-body levels. The low- 
temperature states all have the same intrinsic shape characteristic (the same vortex structure). As we increase 
the temperature, other v phases become accessible. At the point when two bands of different intrinsic character 
cross in energy, we find a phase transition. Similar phenomena are found in nuclear physics, where at higher nuclear 
excitation energies the eigenstates of the system may be of a different intrinsic deformation when compared to states 
of the ground-state band Q . 

We also find generally that the exchange energy decreases as a function of increasing temperature but with differing 
slopes in the different phases. In the region of phase transitions, the ratio of the exchange energy to the direct energy 
decreases from 0.38 (at 11.5 K)to 0.32 (at 13 K) in the v = 6 phase. In the v — 5 phase, this ratio decreases from 
0.32 at T = 13 K to 0.28 at T = 13.8 K, and from 0.28 at T = 13.8 K to 0.26 at T = 14.5 K in the u = A phase. The 
slopes of the decreasing ratio are different across the three phases, with v = 6, 5 being the largest and v — 4 more 
gradual. 

The occupation probabilities of the Fock-Darwin states, ufd, for the temperature conditions of Fig. |l| are shown in 
Fig. ||c. In the low-temperature phase (3.9 K, solid line), the I — state is occupied, while an edge reconstruction has 
occurred for the remaining five electrons. As we increase the temperature to 12 K (dotted line), we see a spreading 
of the occupations in both the low and high angular momentum channels with about 0.8 particles in the I — state. 
The v = 5 phase brings a dramatic decrease of occupation in the I = state and a shift to lower angular momentum 
for the reconstructed edge. This trend continues after the v = 4 transition. 

Signatures of the density transitions that we have seen also appear in the Hartree-Fock occupations uhf as shown 
in Fig. ||d. At low temperatures (3.9 K) the familiar step-function behavior is evident from the figure. At 12 K in 
the v = 6 phase, we see a decrease of occupation to roughly 0.8 in the lowest six Hartree-Fock levels and a spreading 
to higher energy states. As the system undergoes the transition to v = 5, we see only five Hartree-Fock levels 
significantly filled (with uhf > 0.7), and finally in the v = 4 phase, only four Hartree-Fock levels are significantly 
filled. The occupation number-spreading, which is due to thermal excitation of the system, is enhanced significantly 
when the system undergoes a phase transition. 

The phase transitions that we have shown in the preceding discussion have definite observable consequences. In 
Fig. ||a we plot the chemical potential, that is, the separation energy A(N,T) = En(T) — En-i(T) to remove a 
particle from the quantum dot at a given temperature. Since this is an energy difference, A(iV, T) will be influenced 
by transitions within both the N = 6 and N = 5 systems, and we expect changes in slope at the transition points. 
The v = 5 to v = 4 transition in the N = 5 dot occurs at roughly 12 K, causing a sharp rise in A(N = 6). A slope 
change in A (AT = 6) is seen at w 13 K, where the v = 6 to v = 5 transition occurs in the N = 6 system. The decrease 
from 14 — 14.3 K occurs when the N = 6 system makes the transition from v = 5 to v = 4. A final change in slope 
occurs when the N = 5 dot makes the transition from v = 4 to v — 3. For the brief temperature interval when the 
N — 6 and N — 5 dots are in the same v phase, we see a decrease in the chemical potential. 

We observe similar changes in the inverse compressibility, A 2 (N, T) = En+i(T) — 2En(T)+En-i(T). This quantity 
has been measured for quantum dots in low magnetic fields |18| and studied in Hartree-Fock theory for ground-state 
properties [M. In our case, the N = 7,6, and 5 dots participate in the observable. We again notice strong effects as 
one passes through transition points in either of the three systems contributing to the observable. Figure [|b shows 
A2(A = 6) as a function of temperature. Before transitions occur, A2 remains fairly constant. A large decrease 
begins at 12 K, where the N = 5 system undergoes its first transition. Interestingly, A2 increases significantly when 
the N = 6 and N — 5 dots are in the v = 4 phase. 

In order to investigate the sensitivity of our results to the approximations made (namely the Hartree-Fock approxi- 
mation) we extended our studies to include the second-order perturbative correction to the Hartree-Fock energy. This 
correction to the energy is given by 



where we restrict the sums to be below and above the Fermi energy surface, ep. In this expression, the Hartree- 
Fock states are given by a, 6, r, s, with associated single-particle energies e a ,£fc,e r ,e s , and (ab \ v \ rs) are the 
antisymmetrized two-body Hartree-Fock matrix elements of the original two-body interaction of the Hamiltonian. We 
note that the absolute value of AE2 in our calculations is always less than 1% of the total energy, and in fact only 
slightly changes the observable quantities, as is evident from Fig. |^, where the second-order perturbative results are 
also shown for A and A2. 
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In addition to measurements of energy differences, one should be able to experimentally probe the thermal phase 
transitions using far-infrared spectroscopy and X-ray scattering. Far-infrared spectroscopy was used to investigate 
the excitations of InAs quantum dots as a function of the electron number per dot |20|. Grazing incidence X-ray 
scattering was recently used to generate a full structural characterization of quantum dots pl[ , including information 
on the elastic form factor. While these experiments were carried out at very low temperatures, it is conceivable that 
one could study the thermal response of quantum dots using X-ray or far-infrared scattering. A Fourier transform 
of the charge density produced in our calculations gives the elastic form factor that can be used to characterize the 
dot. We show in Fig. || the form factor, | p(q) |, as a function of the momentum transfer vector \ q \= J~q\ + q'y- A 

well-defined minimum is apparent at approximately | q |= 0.7ev _1 in all three phases. This first minimum is related 
to the size of the dot and clearly does not change in the three phases. This is apparent also from a close inspection 
of Fig. || while the internal structure changes significantly as a function of increasing temperature, the dot size does 
not change. The height of the second maximum slightly increases as in phase v = 5 before decreasing in the v = 4 
phase. The position of the second minimum increases as a function of q significantly and would be a distinguishing 
feature in X-ray or far- infrared scattering experiments used to probe the thermal phases of a quantum dot. 

We see in the above calculations that quantum dots exposed to high magnetic fields undergo various phase transitions 
as the temperature increases. Indeed, we have followed the N = 6 system through the v = 3 and 2 phases as well. We 
have also studied these transitions at various magnetic-field strengths, and with differing numbers of electrons in the 
dot up to iV = 20, with the results presented here being characteristic of the general behavior. Our purpose in this 
paper has been to demonstrate that thermal excitations of a quantum dot induce phase transitions within the dot that 
should have observable consequences. Thus, we have shown that temperature may also be used as a control to tune 
the characteristics of quantum dots. Inclusion of the full many-body correlations may slightly dampen some of these 
effects Q , but signatures of the phase transitions should still appear in scattering experiments and measurements of 
the inverse compressibility and chemical potential of the quantum dot. 

This research was sponsored by the Laboratory Directed Research and Development Program of Oak Ridge National 
Laboratory, under Contract No. DE-AC05-00OR22725 managed by UT-Battelle, LLC. by the U.S. Department of 
Energy. 
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FIG. 1. Panels a-d, the charge density; panels e-h, the angular momentum density (in the z-direction) ; panels i-1, the spin 
density. In each panel, —8 nm< x,y < 8 nm. 
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FIG. 2. a) Expectation value of the energy as a function of temperature showing the three phases as discussed in the text, 
b) The specific heat for the lowest-energy configuration of the quantum dot as a function of temperature, c) Occupation of the 
Fock-Darwin states, d) Occupation of the Hartree-Fock states. 
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FIG. 3. a) Chemical potential, A and b) inverse compressibility, A2 for the N = 6 system as a function of the temperature. 
The Hartree-Fock results are shown as a solid line, and the second-order perturbation theory results are given by the dashed 
line. 
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FIG. 4. The form factor in the three phases, v = 6 (11.9 K), v = 5 (13.6 K), and v = 4 (14.3 K). 
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